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Abstract 

We study the stationary points of what is known as the lattice Landau gauge fixing functional 
in one-dimensional compact U(l) lattice gauge theory, or as the Hamiltonian of the one-dimen- 
sional random phase XY model in statistical physics. An analytic solution of all stationary points 
is derived for lattices with an odd number of lattice sites and periodic boundary conditions. In 
the context of lattice gauge theory, these stationary points and their indices are used to compute 
the gauge fixing partition function, making reference in particular to the Neuberger problem. 
Interpreted as stationary points of the one-dimensional XY Hamiltonian, the solutions and their 
Hessian determinants allow us to evaluate a criterion which makes predictions on the existence 
of phase transitions and the corresponding critical energies in the thermodynamic limit. 



1. Introduction 

Stationary points, i. e. points of vanishing gradient of some function, play a prominent role 
in many methods of theoretical physics. Generally speaking, any variational principle uses sta- 
tionarity to single out a certain point in a (function) space. One field of physics where stationary 
points feature prominently is gauge fixing via a gauge fixing functional. Recently, a variety of 
methods summarized under the name of "energy landscape methods" [1] have attracted quite 
some attention, allowing for applications to many-body systems as diverse as metallic clusters, 
biomolecules and their folding transitions, or glass formers undergoing a glass transition. In all 
these examples, the stationary points of some (usually high-dimensional) energy (or potential 
energy, or free energy) function are studied. In most cases such stationary point approaches 
are used in numerical studies. Exact, analytic studies of the stationary points are rarely feasible 
but, as is often the case, can be of great value, as they may prove useful for verifying ideas or 
testing approximation schemes, and they may be instructive towards the understanding of more 
complicated situations. 

In some instances, it may happen that not only the same method of analyzing stationary points 
is used for applications in different fields of physics, but even the same function is analyzed, albeit 
in a very different context. In this article, we want to report such a case, which is an analysis of 
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the stationary points of the function 



N 



(i) 



k=l 



The function F$ is defined on an A^-torus, parametrized by the angular variables 6 = (6\, . . . , %) 
with 6k e (-7T,n]. The random parameters fa e {-jt, n] introduce a quenched disorder, and we 
will be looking for stationary points of with respect to the (^-variables, 



valid for arbitrary fixed values of <p = (fa, . . . , fa/). A central result of this paper is an exact, 
analytic expression for the stationary points of this function for periodic boundary conditions. 
Depending on the problem under investigation, various properties of the stationary points are 
used in the analysis and related to physical quantities. This may be the number of stationary 
points, their respective energy value, the index of a stationary point or its Hessian determinant. 

In physics, the function F$ given in ([T} appears in at least two different contexts: First, it is the 
lattice Landau gauge functional for a compact U(l) lattice gauge theory in one spatial dimension. 
Each of the stationary points, i. e. each solution of the set of equations (fj), corresponds to a gauge 
fixing, and a number of interesting phenomena like the Gribov ambiguity and the Neuberger 
problem are related to the stationary points and their properties. Second, in statistical physics the 
function F$ is known as the Hamiltonian of the classical random-phase XY model in one spatial 
dimension. This model consists of N classical planar spin variables on a one-dimensional lattice 
where each spin is coupled to its two nearest neighbors on the lattice. In recent years, there has 
been significant progress towards establishing a relation between stationary points of classical 
many-body Hamiltonians and the occurrence or absence of phase transitions, and it is in this 
spirit that we will analyze the stationary points of the XY Hamiltonian in this article. Since we 
assume that the average reader will be familiar with at most one of these fields of physics, but 
not both, we will provide elementary introductions in later sections. 

The main concern of the present article is twofold: First, the exact, analytic solution of the 
stationary points of F$ provides a testing ground for ideas, methods, and approximations that 
are used for higher-dimensional or otherwise more complicated models where exact results are 
not available. Second, juxtaposing an interpretation of the stationary points as gauge fixing of 
a compact U(l) lattice gauge theory with an interpretation as special points of a classical spin 
Hamiltonian beautifully illustrates the similarity as well as the differences of stationary point 
approaches in two very different fields of physics. 

The article is structured as follows: In section [2] an exact, analytical expression is derived 
for the stationary points of with periodic boundary conditions. A short introduction to gauge 
fixing in lattice gauge theory, the Gribov ambiguity, and the Neuberger problem is given in 
section [3] In section 13.11 F$ is shown to be the lattice Landau gauge functional of the one- 
dimensional, compact U(l) lattice gauge theory, and the results from section|2]are used to discuss 
the Neuberger problem for this case. Possible strategies to avoid the Neuberger problem and the 
Gribov ambiguity are discussed in section [3~2l In section|4]we turn to the interpretation of F^ as 
the Hamiltonian of a classical random-phase XY spin system. In section |4~T1 some recent results 
on stationary points and their relation to the occurrence of phase transitions are sketched. In 
section l4~2l the random phase XY model is introduced and some thermodynamic properties of the 
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one-dimensional XY model are sketched. In section l4~3l the stationary points of the Hamiltonian 
of this model are put in relation to its thermodynamic behavior. A summary and an outlook on 
open questions and possible future work is given in section[5] 



2. Calculating stationary points 



The stationary points of the function F$ as defined in equation (Q} depend crucially on the 
conditions imposed on the variables 6\ and 6$ at the boundaries. In this section, the stationary 
points of F<f, are calculated for periodic boundary conditions. The stationary points of with 
anti -periodic boundary conditions have been computed previously in yj]. Remarkably, the solu- 
tions found in the two cases differ substantially. For anti-periodic boundary conditions, two fairly 
simple classes of solutions of the set of equations (f2j) were shown in |2[] to exhaustively compre- 
hend the 2 N stationary points of F$, all of which are isolated. For periodic boundary conditions, 
the global rotational symmetry, 6k — > 9k + a for any real a and k = 1, . . .,N, trivially accounts 
for continuous curves of solutions. But even after fixing this global invariance, the solutions for 
the stationary points, as will be shown below, are significantly more variegated than in the case 
of anti-periodic boundary conditions. 

To compute the stationary points of F<f, with periodic boundary conditions, it is convenient to 
transform the set of equations (f2]i to the new angular variables 



(-7T,n] B s k := <p k + i 



6k mod 2n, 



yielding 
together with 



= sin ijt+i — sin s k , k=l,...,N~l, 



= sin si - sin 
Equivalently, we can write these equations in matrix form, 

MS = 0, 



(3) 

(4) 
(5) 

(6) 



with column vector S = (sin s\, . 
matrix of the form 



M = 



. , sin sn) and zero vector 0, where M is a circulant (jV x AO- 
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Like for any circulant matrix, the eigenvalues of M can be readily obtained by a discrete Fourier 
transform of an arbitrary row or column of the matrix 1 3] . One of the eigenvalues is found to be 
zero, corresponding to the invariance of the set of equations (|4]i and ((5]) under the abovementioned 
global phase shift. As a consequence of this invariance, the solutions of (O are not isolated, but 
occur in one-parameter families. Moreover, only N - 1 out of the N equations in © and © are 
linearly independent. To get rid of this trivial invariance, we can set, without loss of generality, 
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On = 0, and eliminate equation (0. Since the inverse function of the sine is multi-valued, we 
obtain multiple solutions for the Sk in ©, 



St + 2nlk = % or st + 2nlk — n - %, (8) 
with 4 being integer valued. Equivalently, we can write 

s k + 2nl k = (-l) 9 *iiv + <Z* 6 (0, 1 }. (9) 

Summing the iV — 1 equations in (O, we obtain 

AT— 1 N-l 

+ 2^ - ffflt) = sjv (10) 

Making use of the definitions S/c = 4>k + - <k for £ = 1, , AT — 1 and % = <Pn + 0\ - #;v> we 

can write 

Af-l JV 

^ j t = % - ft - ^ + ^ <p k = + a>, (ii) 

where we have defined <t> := <Pk- Inserting this expression into ( TTOb and solving for %, one 
obtains 

valid for 1 + Sfj/C-l)* # 0. In case the denominator vanishes, we do not obtain a well defined 
solution for s#. Although we have little doubt that this issue could be resolved, we will simply 
avoid the problem of vanishing denominators by restricting N to odd integers. From the second 
term on the right hand side we can infer that the value of sn does not really depend on the specific 
sequence of h, . . . , In-\, but only on their sum. Without loss of generality, we can therefore write 

Q + 2nl-nZ»:tq k 

S N = rr — : , (13) 

where / is an integer valued parameter. In fact, since we are interested in solutions of % only 
modulo multiples of 2n, it is sufficient to consider 

le (l,..., l+2](-l)* ■ (14) 

k=l 

Summarizing, for N odd, we have obtained an expression for sn, parametrized by the sequence 
(q\, ... ,qN-\) with qic £ {0, 1 } and by the parameter I. To comply with the convention fixed in (O, 
% should then be projected onto the interval {—n, n]. The sequence of variables {s\, . . . , Sn-i) 
is readily obtained by making use of equation (O; here h - for all A: = 1, . . ., N — 1 can be 
used, since we are only interested in solutions for the s^ modulo multiples of 2n. Making use of 
On = 0, it is then straightforward to translate the result to the original variables 0k by inverting 
the definitions ([3). 

Among these solutions, there are two particularly simple classes of stationary points: (a) any 
sequence (si, . .., sn) where s% e {O, O + n} for all k is a stationary point of Ff, and (b) any 
sequence of identical Sk = (O + 2nl)/N with / = 1, . . . , N. 
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The number #(6> s ) of stationary points 8 s that exist for a given number N of lattice sites is 
obtained straightforwardly: Consider all sequences (qi, . . . ,qx-i) consisting of a certain fixed 
number j of 1-entries, and N - j — 1 zero entries. The number of such permutations is (^J 1 ) and, 
according to ( TT4b . the parameter / can take on values ranging from zero to \N - 2 j\. Summing 
over all possible values of j, we obtain the overall number of stationary points, 

W) = |>-2,l( . )' W ^ m . (15) 

where the final expression follows after some tedious but rather straightforward algebra. Not un- 
expectedly Jj], asymptotically for large system sizes N, applying Stirling's formula to ( TT3b shows 
that the number of stationary points increases exponentially with N. The number of minima of 
F$ can be obtained by reverting to the following result. 

Theorem 1. Consider the function F$ as defined in ([T]) and its stationary points as derived 
above, each of which is parametrized by a sequence {q\, . . . , q^) with q^ e {0, 1 } and by I € 
{1, . . ., 1 + 2jfc=i ( — l) ?t }- A stationary point of 'F^ is a minimum if and only if 

(a) qt — Ofor all k — 1, ... ,N — 1, and 

(b) cos sn > 0, with s^ as defined in (113 b . 

The proof of this result is similar to a related statement shown in the Appendix of J^t], and 
we refer the reader to this article for more information on the strategy of proof. Note that, via the 
O-dependence of $n, the number of minima of F<f, depends in general on the quenched disorder 
4> = (0i, ... , fa). Interestingly, one can read off from the above theorem that the number of 
minima is bounded from above by N. The number of minima increases therefore at most linearly 
with N, but not exponentially, in contrast to the general expectation based on an assumption of 
independent subsystems yj. 

Evaluating the function ([1]) at a stationary point, we obtain 

N 

F m = N -cos s N J](-l) q K (16) 

k=\ 

The Hessian = CHij) of F^ is a tridiagonal matrix with elements 
, d 2 FJ0) 

7 = 86 88- = 6 i,j( coss i-i +coss i)- 6 i-hj coss i-i- 6 MJ coss ^ 'J = 1, - . .,iV— 1. (17) 

Note that H is an (N - 1) x (N - l)-matrix, accounting for the fact that we have removed one 
linearly dependent equation out of the equations © and ((5]). The Hessian determinant can be 
computed explicitly, yielding 

N N 



i=\ k=\ 

Evaluating *H at the stationary points s of F$, we find 

detW) = (cos^) JV - 1 ri ( " 1)9 ' ' (19) 

\ £-1 )\ j =l I 



3. Lattice Landau gauge fixing 



In this section, a brief introduction to gauge fixing on the lattice is given. In section 13.11 
the function defined in (Q]) will be identified as the Landau gauge fixing functional of a 
one-dimensional compact U(l) lattice gauge theory, and this connection allows us to use the 
stationary points of computed in section|2]for studying the gauge fixing partition function. 

A gauge field theory can be successfully studied nonperturbatively by discretizing the Eu- 
clidean space-time and putting the gauge and matter fields of the theory on a (d+ l)-dimensional 
space-time grid [6]. A lattice grid consists of lattice sites, as well as links connecting any two ad- 
jacent sites. The gauge fields are usually defined through link variables E/; >jU e G where i denotes 
the site index, fi is a directional index, and G is the gauge group of the theory. By U — {t/,^} we 
define the set of link variables for all i and ft. These link variables are related to the gauge fields 
of the continuum theory by the standard Wilson discretization |6|]. 

Gauge fixing is a procedure to get rid of the redundant degrees of freedom of the field vari- 
ables in gauge theories. Technically, gauge fixing is implemented by requiring the gauge fields to 
satisfy a certain constraint, like for example one of the well-known gauge conditions in classical 
electrodynamics. In a continuum gauge theory, gauge fixing is indispensable in order to assure 
convergence of certain integrals when computing expectation values of observables. In lattice 
gauge theory, the spatial discretization renders the theory manifestly gauge invariant: when com- 
puting expectation values of an observable O with respect to the (gauge-invariant) Wilson action 
S, 

= jn^du^xpi-smwiu] 

fUi,,dU^cx P (-S[U]) 

the occurring integrals are well-defined, and fixing a gauge is therefore not required. Still, fix- 
ing a gauge on the lattice can be beneficial, for example when comparing with gauge dependent 
results from the corresponding continuum theory. For example, in continuum gauge theories, 
gauge-dependent quantities are at the basis of certain confinement scenarios such as Kugo-Ojima 
and Gribov-Zwanziger J3,Htl, or Dyson-Schwinger equation studies for the quantum chromody- 
namics hadron phenomenology ^ . 

In its modern formulation, a gauge fixing condition is usually specified in terms of the minima 
or, more generally, the stationary points of some functional I2l ll0l4l2ll . On the lattice, a popular 
choice is the lattice Landau gauge functional [ 13], 

Fu(g) = 2 t 1 - Re T*glUi*gi#) . ( 21 ) 

where g = (gi,g2, ■ ■ •) with gj E G are gauge transformations, and the summation in (fJTJ runs 
over all sites i and all link directions fi on the lattice. The gauge-fixing conditions are then defined 
as the vanishing first derivatives of Fu with respect to the gauge transformations g, 

/Qr) = C/i(g),/2(g),...) = with fi {g):= d I^-. (22) 

dgi 

The functional (I2TI) is chosen such that the gauge-fixing equations (l22l) reproduce the usual con- 
tinuum Landau gauge in the naive continuum limit. 

When computing expectation values of some observable O, gauge fixing is implemented by 
constraining the integrations to those values of the link variables U lp which satisfy the gauge 
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fixing equations d22l) . Technically this can be done by a procedure named after Faddeev and 
Popov 11411 . resulting in the expression 

= [ 11^ dU^A(g)6(f(g)) exp(-S [U])Q[U]/Z GF 
J n^dU^AigWfig)) exp(-S [t/])/Z GF 

for the expectation value of a gauge invariant observable 0[U], The so-called Faddeev-Popov 
determinant A = det "H in this formula is the determinant of the Hessian matrix *H of Fjj with 
elements 

%) = (24) 

evaluated at the stationary points g s of Fjj. The gauge fixing partition function occurring in the 
expectation value (1231 ) can be written in the form 

z GF = 2 s § nA ^ s )' < 25 ) 



where the summation is over all stationary points of Fy 115 ] . 

In the 1970s, Gribov noticed that gauge fixing can in general not remove the ambiguity in 
the field variables completely [8]. Instead, even after fixing a local gauge, there exist gauge- 
equivalent field configurations. In our setting, these so-called Gribov copies manifest themselves 
in multiple stationary points of the gauge fixing functional OTI) . Despite the Gribov ambiguity, 
the Faddeev-Popov procedure for deriving the gauge-fixed expectation value works fine in the 
continuum, but a serious problem arises for compact U(l) and SU(N C ) gauge theories on the 
lattice: Neuberger observed in the 1980s that, in these cases, the Gribov ambiguity leads to a 
vanishing lattice gauge fixing partition function Zgf 116, 17I1 . and hence to vanishing denomi- 



nators in the gauge-fixed expectation value d23l ). From this so-called Neuberger problem, one 
has to conclude that the standard lattice Landau gauge OTT) is not a suitable gauge fixing condi- 
tion for compact U(l) and SU(N C ) gauge theories on the lattice. This failure can be viewed as 
originating from a deeper cause: for compact U(l) or SU(iV e ) gauge groups, the standard lattice 
Landau gauge fixing procedure cannot restore a Becchi-Rouet-Stora-Tyutin (BRST) symmetry 
on the lattice. 

In fact, Neuberger first observed a vanishing Zgf when he tried to establish BRST symmetry 
on the lattice. Later, Schaden interpreted the problem in terms of Morse theory and showed 



in [15] that, by virtue of the Poincare-Hopf theorem, one can write the gauge fixing partition 
function as 

Z GP = X (G N ) = \ X {G)f , (26) 

where x(G) denotes the Euler characteristic of the group manifold G. The Euler characteristic 
is a topological invariant, and for both, U(l) and SU(N C ) group manifolds, we have x(G) = 0, 
and Zgf therefore vanishes. Making use of such a topological interpretation, it was argued in 
that the Neuberger problem for an SU(N C ) lattice gauge theory actually lies in the compact U(l) 
subgroup of S\J(N C ). Therefore, finding a solution to the U(l) Neuberger problem is expected 
to have an even bigger impact, as it should allow to evade the problem also for SU(N C ) gauge 
groups. 

In the next section, we will use the stationary points computed in section [2] to explicitly 
demonstrate the Neuberger problem for a compact U(l) gauge theory on a one-dimensional lat- 
tice. 
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3.1. One -dimensional compact U(l) lattice gauge theory 

Topologically, the group manifold of a compact U(l) gauge group is a circle. We parametrize 
this circle by means of angular variables, the link angles (f>i tU e (-n, n] and the gauge trans- 
formation angles € {-n, jt]. In terms of these angles, we can write the link variables as 
t/, ;/J = exp(i0 1AI ) and the gauge transformations as g,- = exp(i6',). Inserting these expressions 
into equation ( f2TT > and specializing to a one-dimensional lattice, the lattice Landau gauge fixing 
functional is identical to the function F<p defined in ([TJ. This allows us to use the stationary points 
and Hessian determinants obtained in section|2]for computing the gauge fixing partition function 
(l25T l of a one-dimensional compact U(l) lattice gauge theory. 

For even N, the computation of Zqf is particularly simple and does not even require explicit 
solutions for the stationary points: For every stationary point s = (si , . . . , s/y) satisfying the gauge 
fixing equations © and (0, there is a gauge transformed Gribov copy s' — (s\ + n, . . . , Sn + n), 
likewise satisfying the gauge fixing equations. From equation (TT~8b we can infer that the Faddeev- 
Popov determinants (or Hessian determinants) at s and s' have opposite signs, 

N N N N 

detH(s') = ^ Y\ cos( ^ + n ) = _ X Ft cos Sk = - Aet<H ^)- ( 2? ) 

7=1 k=l i=i k=\ 

k*j k*j 

This symmetry allows us to easily read off that the gauge fixing partition function 

Z GF = ^sgnA(0 s ) = O (28) 

vanishes for the one-dimensional compact U(l) lattice gauge theory with an even number N of 
lattice sites. 

For odd N, no symmetry of the type (l27t exist. In order to evaluate the gauge fixing partition 
function (l25T l in this case, we make use of the Hessian determinant evaluated at the stationary 
points ([T9T l, as computed in section|2] Since (cos S/y)^ -1 ' is positive for odd N, the sign of ([T9T l 
is determined by the remaining terms 

n (_i) *i- (29) 

This expression evidently does not depend on the parameter / labelling a particular stationary 
point, nor on the specific sequence of binary parameters e {0, 1 }, but only on the number 



-2 



N-l 

q k (30) 

k=l 



of "l"s in such a sequence. Writing (1291 1 in terms of n, we obtain (N - 2n)(-l)", and the sign of 
this expression determines the sign of the Hessian determinant at a given stationary point. Since 
there are sequences (qi, . . . , qN-\) with a given n, times \N - 2n\ different parameters /, we 
obtain 



N-l I j, y. N-l , y. 

Z GF = J]sgnA(0 s ) = 2 ~ W-2«|sgn(^V-2«)(-l)" = ^ ~ \{N - 2n)(-l)" = 0. 



(31) 



This result illustrates the occurrence of the Neuberger problem, i.e. of a vanishing gauge fixing 
partition function, for a one-dimensional compact U(l) lattice gauge theory. Since this leads 
to vanishing denominators in expectation values, and therefore to meaningless expressions, we 
must conclude that the chosen gauge fixing procedure is not adequate for this gauge theory. 

3.2. Modified lattice Landau gauges 

In order to avoid the Neuberger problem, we might want to modify the gauge fixing condition 
(f2]i in some way. Apart from avoiding the Neuberger problem, there are further properties a mod- 
ified gauge fixing condition should satisfy. In particular, the computational complexity should 
be such that a numerical implementation is feasible in practice, and the gauge fixing partition 
function Zgf should be independent of the orbit 0. A few possible modifications of the lattice 
Landau gauge are discussed in this section. 



Minimal lattice Landau gauge. We had seen in sectionQJJthat Zgf vanishes due to cancellations 
of Hessian determinants with opposite signs. An obvious strategy to avoid such cancellations 
consists in strengthening the gauge fixing constraint, for example by requiring not only a sta- 
tionary point of the gauge fixing functional ([T), but a local minimum |8J]. Apart from avoiding 
the Neuberger problem, this so-called minimal lattice Landau gauge is also numerically advanta- 
geous, as minima are easier to track down than general stationary points. However, as observed 
in theorem[T]and the subsequent discussion, the number of minima depends on the gauge orbit <p. 
Hence, also the gauge fixing partition function Zgf is orbit-dependent, and this is an undesired 
effect. 

Absolute lattice Landau gauge. The unwanted dependence of Zgf on the gauge orbit can 
generically be eliminated by further strengthening the gauge fixing constraint to the set of global 
minima of the gauge fixing functional (Q]), and this choice is called absolute lattice Landau gauge. 
From ( TToT l and theorem[T]we find that the absolute minima of F^, are given by the stationary points 
with qk — and a parameter I chosen such that cos[(<E> + 2nl)/N] is maximized. 

Although the absolute lattice Landau gauge resolves both the Neuberger problem and, for 
generic values of <t>, the problem of orbit-dependent Zgf, it suffers from another drawback: In 
higher dimensions, where no analytical results are available, finding the global minimum of the 
gauge fixing functional is a computationally very hard problem, rendering an application of this 
gauge fixing condition practically impossible. 

Stereographically modified lattice Landau gauge. Another attempt to overcome the Neuberger 
problem, based on a stereographic projection of the U(l) group manifold, has been proposed 
and studied recently iB 1 1 OL 121. The stereographically modified lattice Landau gauge fixing 



functional for the compact U(l) case in one dimension can be written as 

PM = -2 2 In ( l+c°»(fr + ft+.-ft> ) . (32) 

For this function one can show that the number of stationary points is independent of the orbit 
<p and exponentially (in N) smaller compared to the standard lattice Landau gauge [18]. More- 
over, after removing the global gauge freedom, the Hessian determinant is generically positive 
definite, cancellation of signs of determinants cannot occur, and the Neuberger problem is there- 



fore absent I12L 11911 . This latter property also generalizes to gauge groups SU(JV C ) and to higher 
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dimensional lattices. It remains, however, an open question whether the desired property of 
orbit-independence of Zgf generalizes to higher-dimensional versions of the functional (|32| |. In 
summary, the stereographically modified lattice Landau gauge satisfactorily solves the Neuberger 
problem on a one-dimensional lattice, but it remains to be seen whether this is true also in higher 
dimensions. 



4. Random-phase XY model 

In this section, some recent results are reviewed that relate stationary points of the Hamil- 
tonian function of a classical many-body system to the occurrence of a phase transition in the 
thermodynamic limit of large system size. The function defined in (Q]) is identified as the 
Hamiltonian of the random phase XY model in one spatial dimension. This connection allows 
us to use the stationary points of F^ and their Hessian determinants computed in section[2]for an 
analysis of certain thermodynamic properties of the model. 



4.1. Stationary points and phase transitions 

As is long known, the stationary points of a classical Hamiltonian function can be employed 
to calculate or estimate certain physical quantities of interest. Famous examples include transi- 
tion state theory 112011 or Kramers's reaction rate theory for the thermally activated escape from 
metastable states 12111 . where the barrier height, i. e. the potential at a certain stationary point of 
the potential energy function, plays an essential role. More recently, the noise-free escape from 
quasi-stationary metastable states, whose lifetimes diverge with the system size, has been related 
to the presence of stationary points of marginal stability Ir22ll . Apart from studies of dynami- 
cal properties, stationary points have also been extensively used for estimating thermodynamical 



properties by means of the superposition approach M23I1 



Dynamical properties like the aforementioned ones are, as one might expect, not unrelated 



to the statistical physical behaviour of a system. Accordingly, as worked out beautifully in 112411 . 
properties of stationary point^] reflect in dynamical and statistical physical quantities simultane- 
ously. This observation sparked quite some research activity, reviewed in 12511 . with the aim of 
relating equilibrium phase transitions to stationary points and their indices. Among other things, 
it was observed that stationary points of the potential energy function generate nonanaiyticities 
of the microcanonical entropy like the ones reported in B26I1 . Subsequently, it was noticed in 
1 27 - 29| that the Hessian determinant of the Hamiltonian, evaluated at the stationary points, adds 



a crucial piece of information that helps to discriminate whether or not a phase transition occurs. 
Omitting the technical details, the essence of the criterion on the Hessian determinant can be 
captured as follows 1 30] . 

Criterion. Let 

1 N 

H{p,q)=-Y J p 2 k + V( gi ,...,q N ) (33) 
k=i 

be the Hamiltonian of a system with N degrees of freedom, possessing a number of stationary 
points that grows at most exponentially with N. Here, p — (p\,,.. ,Pn) and q — (q\, . . . ,q^) 



1 The authors of [ 24] discuss topology changes of constant-energy manifolds in phase space. Via Morse theory, such 
topology changes can be related to the stationary points of a sufficiently smooth Hamiltonian. 
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denote, respectively, the vectors of momenta and positions, and V is the potential energy. In the 
thermodynamic limit N — > oo, a phase transition can occur at some critical energy per degree of 
freedom e c only if the following two conditions are met: 

1. There exists a sequence {q^}^ =No of stationary points ofV such that 

Vc:= j™ — (34) 

converges and equals (v)(e c ) (i.e. the ensemble expectation value ofv — V/N at the energy 
e c ). 

2. The asymptotic behaviour of the Hessian matrix 'HofV, evaluated at the critical points qf 
contained in that sequence, is such that 

limldet 7/(^)1 1/A? = 0. (35) 

Note that this criterion is necessary, but not sufficient, for a phase transition to occur: finding 
a sequence of stationary points with the behaviour specified above does not guarantee a transi- 
tion to take place at the corresponding critical energy. However, as model calculations show, 
the criterion usually appears to single out precisely the correct transition energies 11281 12911 . Im- 
portantly for the application of the above criterion, knowledge of a suitably chosen subset of all 
the stationary points of V may be sufficient. This matter of fact was pointed out by Nardini and 
Casetti [30], and suitably constructed sequences of stationary points were used to single out the 
phase transition and determine its critical energy for a model of gravitating masses. 

Comparing the above criterion to other analytic tools in the statistical physics of phase tran- 
sitions, a remarkable property is its locality in configuration space. In contrast to, say, the calcu- 
lation of a partition function, no averaging over a large, high-dimensional manifold is necessary. 
Instead, only the local properties of a sequence of stationary points needs to be analyzed. Of 
course, finding an appropriate sequence of stationary points can be equally hard or impossible, 
but in certain instances such a local approach may prove beneficial. 

4.2. Random phase XY model 

The study of spin models has been very fruitful towards the understanding of phase tran- 
sitions. First, such models are typically simpler than continuum systems, and second they are 
successfully employed for the modelling of condensed matter phenomena like ferromagnetism 
and others. An exceptionally popular class of such models are O(n) vector models, consisting 
of n-dimensional unit vectors crk placed on the sites of a lattice. The interaction between the 
classical spin vectors is described by a Hamiltonian of the form 

H{o-) = jYj^-<rk-o-i) (36) 

(k,V> 

with coupling constant /. The angular brackets (■, ■) denote a summation over all pairs of nearest- 
neighbouring sites on the lattice. O(n) vector models include as special cases (a) the Ising model, 
n = 1, (b) the XY model, n = 2, and (c) the Heisenberg model, n — 3. For n — 2, we can 
parametrize the two-dimensional unit vectors cr^ by an angular variable 6%, obtaining the Hamil- 
tonian of the classical XY model, 

H(0) = jY J [l-cos(6 k -Oi)l (37) 
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A number of disordered variants of 0(n) vector models can be found in the literature, and 
their study has improved the understanding of glassy behaviour, ageing, and related phenomena. 
One such disordered variant, the random phase XY model, is characterized by the Hamiltonian 



H <p (0)=J2Jl-cos(cf> l + 6 k -6i)l (38) 

{k,D 

where the (pi are angular variables, drawn randomly according to some probability distribution. 
On a one-dimensional lattice and with coupling / = 1, this Hamiltonian takes on the form of F,p 
as defined in (Q}. 

On a lattice of spatial dimension d — 2 or higher, the random phases 0j cause frustration, 
resulting in a number of interesting physical properties. In one dimension, however, frustration 
does not occur, and the random phases do not make much of a difference compared to the ordered 
model (f3Tb . We will therefore mainly focus on the one-dimensional XY model in the absence of 
disorder. The Hamiltonian of this model is 

N 

H{0) = F o (0) = - cos(ft +I - fa], (39) 
*=i 

i. e. as defined in (Q} with 0* = for all k = 1, . . . , N. An exact solution for the canonical 
partition function (or, equivalently, for the canonical free energy /) as a function of the inverse 



temperature f3 has been reported in [31]. In the thermodynamic limit N — » oo and translated into 
our notation, the solution is given by 

/0) = l-±ln/ o <0), (40) 

where Io denotes the modified Bessel function of the first kind. For our purposes, the entropy 
s as a function of the energy e is more convenient, and can be obtained from / by means of a 
Legendre transform, 

s(e)=fi{e-Y)+\nime)), (41) 

where /?(e) is the solution of 

hip) = (1 - e)/o0§). (42) 

The energy e can take on values from the interval [0,2]. The graph of s is shown in figure Q] 
together with a plot of the temperature Tie) = [s'(e)]~' . 

Although this model has no phase transition, it will prove instructive as what regards the 
study of stationary points of many-body systems in the spirit of section |4~T1 

4.3. Stationary points analysis of the one-dimensional XY model 

To study the one-dimensional XY model in the spirit of the criterion stated in section |4~T1 we 
need to compute, for suitable sequences of stationary points s , the (potential) energy e^iff*) = 
Foi(P)/N together with the scaled Hessian determinant, 

£V(0 S ) := |detW(fl ,, )| 1/(Af_1) . (43) 

Inserting the stationary points specified by (0 with sn given by ( TTSl into the Hamiltonian ( 1391 , 
we obtain 

j , N-l 

e N iO s ) = 1 - -(l + 2(-l)«Jcos s N . (44) 
12 
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Figure 1: Graph of the entropy * and the temperature T, both as functions of the energy e, for the one-dimensional XY 
model. The range of accessible energies is [0, 2], where energies < e < 1 correspond to positive temperatures, energies 
1 < e < 2 to negative temperatures. 



For the scaled Hessian determinant, we find 



£V(f/) = |coss w |g](-l)* 



l/Ctf-i) 



(45) 



by making use of the expression ( fl9] l for the Hessian determinant. 

To apply the criterion of section |4~T1 we need to choose appropriate sequences {# s }^ =iVo °f 
stationary points such that the corresponding sequence of energies e^(f? s ) converges. This is 
most conveniently done by fixing values of and of the parameter 



[-l,l]3/ ? := I(l +£(-!)•). 



For any sequenc^l of stationary points with fixed s# and l q , we can write 

<?Ar(# S ) = I ~lq COS SjV 

and 

D N (^) = \coss N \\Nl q -l\ 1/(N - 1) . 
In the thermodynamic limit N — > oo, the latter simplifies to 

£) := lim £>,v(6> s ) = |cos jjvl, 

;V— >co 



(46) 

(47) 
(48) 

(49) 



becoming independent of 4. Note that, by definition, is fixed for the sequence of station- 
ary points we are considering and is therefore not affected by the limiting procedure in d49b . 
The possible values of s^, as determined by ( TPl ). become dense on the interval {—n, n] in the 
thermodynamic limit. Combining equations ( l47T i and j49l , we obtain 



1 - e N (0 s ) 



la 



(50) 



This implies that, for all stationary points f?, the pairs {e^,!)) are located within the shaded 
triangle plotted in figure|2] and in fact they fill the triangle densely in the limit — ♦ oo. 



2 To construct a sequence of stationary points with a fixed value of sy, it is in general necessary to restrict the sequence 
to some infinite subset of system sizes N for which solutions with the chosen value of sn exist. 
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0.0 0.5 1.0 1.5 2.0 

e 

Figure 2: In the plane of energy per degree of freedom e and scaled Hessian determinant D, the region is shown for 
which stationary points 6 s with the corresponding values (e, D) exist. In the thermodynamic limit N — > oa, the stationary 
points fill the shaded region densely. 



Interpreting these results with regard to the criterion stated in section I4TTI we first have to 
note that, in contrast to ( 1331 ). the Hamiltonian of the one-dimensional XY model has no kinetic 
energy term \ Yll=\ p\- Luckily, this is of minor importance and the criterion still applies, with 
the only simplification that critical energies e c and critical potential energies v c coincide (see 13211 
for a related discussion). 

From the fact that T> > for all e + 1 we can conclude that the one-dimensional XY model 
has no phase transition for these values of the energy. Only for e = 1, the criterion cannot exclude 
a phase transition. Of course, we know from the thermodynamic solution that no such transition 
exists. As pointed out in section |4~T1 the criterion on the Hessian determinant allows to exclude 
a phase transition when T> is bounded away from zero, but a vanishing T> does not necessary 
imply the occurrence of a phase transition at the corresponding energy e0 Moreover, note that 
e — 1 is a very peculiar energy value of the one-dimensional XY model: This energy corresponds 
to temperatures T = +oo, as is evident from the right plot in figureQ] So even if the entropy s(e) 
had shown a nonanalyticity at e = 1, the corresponding phase transition would have occurred at 
infinite temperature. 

Although we could not expect much news on the well-known thermodynamics of the one- 
dimensional XY model, there are a number of other interesting observations to be made: The plot 
in figure[2]shows a remarkable difference compared to earlier work in which phase transitions had 
been analyzed via the Hessian determinant at stationary points: For the mean-field type models 
studied in ll28j|29i[33ll . a unique value of D was found for each value of the energy e. In contrast, 
for the one-dimensional XY model, we observe ©-values densely covering an interval. This 
property, which we assume to be generic, has remarkable consequences for the application of 
the Hessian determinant criterion: It is, apparently, not sufficient to study any sequence {# S }^ =A , 
for which the sequence of energies e^(f? s ) converges. Instead, in order to single out energies for 
which T> can become zero, a suitable sequence of stationary points has to be employed. To get 
an idea what sequences are the relevant ones, the following considerations are of help. 

In section [2] we had remarked that, among the many stationary points of the function F^,, 
there exist two classes of particularly simple ones: (a) The solutions (si, . . . , s^) with Sk e {0, n) 
for all k. These solutions, by virtue of equation (1491 1. all yield !D = 1, and they densely populate 
the upper boundary of the triangle in figure [2] (b) A second class of simple solutions consists of 
identical values = 2nl/N for all lc, where I is some integer between 1 and Af. These solutions 
all correspond to — for all k — \,...,N—\, and therefore l q — 1 . By means of equation ( 



3 A vanishing D in the absence of a phase transition had been observed earlier for a mean-field model of two particle 
species 13311 . 
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the solutions are found to be located on the two equal edges of the isosceles triangle in figure [2 
and they provide lower bounds on D for all accessible energies e e [0, 2]. 

Hence, for the purpose of singling out those energies for which the rescaled Hessian determi- 
nant T> can be vanishing, this second class of solutions with Sk = 2t:1/N for all k would have been 
sufficient. What is more, these symmetric solutions are particularly easy to analyze. Maybe this 
points towards a deeper reason behind the remarkable success of Nardini and Casetti's analysis 
of a model of gravitating masses on a ring, using only a subset of symmetric stationary points 
& 

Another interesting observation arises when comparing results for the stationary points of F$ 
obtained for different boundary conditions. For anti-periodic boundary conditions, s^ + \ = —s\ 
and so = ~ s n, a calculation of the stationary points has been reported in ]2|,[l2|]. In this case, all 
2 N stationary points are of type (a), i.e. (s\, . . ., %) with Sk e {0,7r). Plotting Das a function of 
the energy e, only the horizontal line D = 1 in figure[2]is obtained. According to the criterion of 
section 14.11 this result rules out the existence of a phase transition in the thermodynamic limit. 
Since thermodynamic quantities of short-range interacting systems are known to be independent 
of the boundary conditions, we can infer the absence of a phase transition for the one-dimensional 
nearest-neighbour XY model with any boundary conditions from the stationary points in the case 
of anti-periodic boundary conditions. Again, this observation might prove useful when studying 
other models: Since the boundary conditions strongly affect the stationary points, a suitable 
choice of the boundary conditions might simplify the computation in some cases, and/or allow 
to establish the criterion of section l4~T1 in others. 

Let us finally remark that stationary points of the one-dimensional XY model and their rela- 
tion to phase transitions had been analyzed previously by Casetti et al. in section 3 of (5J. There 
are, however, a number of important differences to our work. First, these authors employ a dif- 
ferent strategy to deal with the global invariance 6k — > 6k + a. Instead of fixing one variable 
0n = as in our analysis, they introduce a symmetry-breaking external magnetic field in the 
Hamiltonian. Second, the study in [5] focusses on stationary points and their indices, but does 
not consider the Hessian determinant at these points. This is not much of a surprise since, at 
that time, the importance of the Hessian determinant at the stationary points had not yet been 
realized. Third, Casetti et al. consider those special stationary points where Sk e {0, n} to be 
exhaustive, missing out the many other ones. Considering this subset of stationary points only, 
the vanishing rescaled Hessian determinant D at e = 1 would have gone unnoticed. 

5. Summary and Outlook 

We have studied the function F^,{6) — Zf =1 [l - cos(<pk + 6k+i - 6k)] with periodic boundary 
conditions, obtaining all its stationary points for any odd N. We have employed this result to 
illustrate the relevance and usefulness of stationary points in different branches of physics. 

When interpreting F$ as the lattice Landau gauge fixing functional in one-dimensional com- 
pact U(l) lattice gauge theory, we used the exact solution for the stationary points 6 s (called 
Gribov copies in this context) and their Hessian determinants (or Faddeev-Popov determinants) 
A(8 S ) to show that the gauge fixing partition Zgf = E s g n A(0 S ) vanishes. The Neuberger problem 
of vanishing Zgf results in ill-defined expectation values, pointing towards a deeper problem: the 
standard lattice Landau gauge fixing fails to establish the BRST symmetry on the lattice, and is 
therefore not adequate for the one-dimensional compact U(l) lattice gauge theory. Also based 
on the exact expression for the stationary points, a number of strategies of how to avoid the 
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Neuberger problem and the Gribov ambiguity were discussed. The hope is that a thorough un- 
derstanding of the one-dimensional case is helpful for devising strategies of how to avoid these 
problems also for gauge theories on higher dimensional lattices. 

Interpreting F$ as the Hamiltonian of the one-dimensional random phase XY model in clas- 
sical statistical physics, the solutions for the stationary points and their Hessian determinants 
allowed us to evaluate a criterion which makes predictions on the existence of phase transitions 
and their critical energies in the thermodynamic limit. In particular, unlike the long-range inter- 
acting models analyzed earlier, the computation provides the first example of a model where, for 
a given value of the energy e, the thermodynamic limit D of the Hessian determinant densely 
covers a range of values. Furthermore, we observed that the boundary conditions have a drastic 
effect on the stationary points, and therefore on the values of D. These findings, we believe, can 
be helpful when studying more complicated models where a complete solution for the stationary 
points is not feasible. 

For lattices of dimension two or larger, the XY model is known to undergo a phase transition. 
The XY model on a two-dimensional square lattice therefore appears to be an ideal choice for 
generalizing the methods and results of the present article. In particular, solutions of the form 
(0 straightforwardly generalize to higher-dimensional lattices. However, for dimensions larger 
than one they cannot be expected to be the only solutions of the stationary point equations, but 
many others might exist. Apart from such an analytical approach, one of us is currently working 
on numerical approaches for computing the stationary points of XY models in dimensions larger 



than one 11341 135H . In order to check whether all stationary points have been found in such a 
computation, one can use the Neuberger zero as a necessary, though not sufficient, condition. 

Apart from the interpretation of F$ as a Hamiltonian function, its stationary points can also be 
used to discuss non-Hamiltonian dynamical systems: The stationary points of Fq can be shown to 
be the fixed points of the one-dimensional Kuramoto model with nearest-neighbour interactions, 

de k K r - - i 

— = a) k + — |sin(6k+i - k ) - sm{e k - 9 k -i)\ , (51) 

provided the frequencies cj k — a> are identical for all k = l,...,N. Under the latter condition, 
one can transform the equations of motion to new variables 8 k =6 k — ajt, yielding 

~i7 = - [sin(0 t+I - k ) - sm(6 k - 6t-i )] , (52) 
at 2 

where K e 1R is some coupling constant. At a stationary point s of F o, the right-hand side 
of d52b vanishes, resulting in a fixed point of the dynamics. The fixed point is stable if s is a 
minimum of Fq as characterized in TheoremQ] Although the one-dimensional nearest-neighbour 
Kuramoto model is known to show no synchronization transition ll36Tl . its fixed points may be 
useful for studying synchronization-related dynamical features of finite systems 13711 . 

We conclude this outlook by mentioning another example of a lattice gauge fixing functional 



that can be interpreted as a physical system: It was observed in 1 38] that the lattice gauge fixing 
functional of the linear covariant gauge can be interpreted as the Hamiltonian of a spin-glass 
model in a random external magnetic field. For a compact U(l) gauge group, the stationary 
points of this gauge fixing functional also correspond to the fixed points of a Kuramoto model 
with nearest neighbour interactions and frequencies a> k distributed according to some probability 
distribution. For other compact Lie groups, similar relations are expected to exist between the 
lattice gauge fixing functionals of the linear covariant gauge and non-Abelian generalizations 
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of the Kuramoto model [39]. Similar to the results reported in the present article, a stationary 
point analysis may lead to interesting cross-connections linking the lattice gauge and dynamical 
system interpretations of these models. 
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